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1. Introduction 

The one-center location problem in the plane, also 
called the one facility min-max location problem, may 
be stated as follows: given m points/?, e IE?, i = 1, 2, ..., 
m, and some distance function d{x,y) for x, ye R^, the 
problem is to determine the location of a point x e R^ 
that minimizes the maximum distance dip^, x) over / = 
1, ..., m. The problem is denoted by PI and written as 
follows: 



PI : min max d{p., x). 

xelR^ i=l,...m 



The equivalent, constrained version of PI is written as 
follows: 

PI: min z 

s.t. z > d{pi, x) i = 1, ..., m. 

Problem PI was first reported by Sylvester [6] using 
Euclidean distance. A substantial number of papers 



have appeared on the one-center problem assuming a 
variety of distances. Articles and books with reviews of 
the literature include Heam and Vijay, [4], Drezner [1], 
Drezner and Hamacher [2]. 

A related problem is the one-median, or total-cost, 
location problem in which the maximum operator is 
replaced by the summation operator. Both the center 
and the median problems have been studied extensive- 
ly for Euclidean distance, for /^ distances for 1 <jt? < oo, 
and on networks [1,2]. Witzgall [9] considered the 
median problem for polyhedral norms and noted that 
the problem could be formulated as a linear program- 
ming problem. He also noted that "the linear program 
had special properties that should be exploited for an 
efficient solution", and that "more research in this area 
is indicated." 

Ward and Wendell [8] considered both the one-medi- 
an and the one-center problems using block distance, 
which is the special case of polyhedral distances in 
which the polytope is symmetric. For both the one- 
median and the one-center problems they reported two 
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linear programming formulations based on characteri- 
zations of block distance in terms of fundamental direc- 
tions and polar directions. 

This paper considers the two linear programming 
formulations of the one-center problem with block dis- 
tances, as presented by Ward and Wendell. The equiva- 
lence of these two formulations follows from the equiv- 
alence of the block distance representations. The dual 
simplex algorithm is applied to the linear programming 
formulation based on polar directions of the block dis- 
tance, and a geometric interpretation is presented. This 
interpretation is applied to the Euclidean distance one- 
center problem and provides an alternative update pro- 
cedure for the dual algorithm. 

This paper actually considers a generalization of P 1 , 
the weighted one-center problem, in which there is a 
positive weight Wi associated with each point/?;, / = 1, 
..., m. The problem is denoted by P2 and the con- 
strained version is stated as follows: 

P2: min z 

s,t z > Wfi(pi, x) i = 1, ..., m. 



2. Block Distance 

Block distance is a special case of general norms and 
were introduced to location problems by Witzgall [9] 
and by Ward and Wendell [7,8]. Block distance is 
defined in the plane with respect to a symmetric poly- 
tope as its unit ball, denoted by B. The polytope B is 
assumed to have 2p distinct extreme points, for some 
integer/? > 2. The vectors corresponding to the extreme 
points are cMcd fundamental directions, and are denot- 
ed by ^1, bj, ..., bj^ where b^+k = -bk for k= \, ..., /?. 
Assume that the fundamental directions are ordered 
counter clockwise, and for notational convenience, let 
bip^k = bj,for k = 1, ...,/?. Figure 1 shows an example 
with p = 4 fiindamental directions and the correspon- 
ding unit ball. 

The block distance between the points x^ and x^ with 
respect to a given set of 2/> fiindamental directions bi, 
..., b2p, is denoted d^ix^, x^) and is defined to be the 
objective fiinction value of the following linear pro- 
gramming problem, denoted by LPD: 

LPD: d^(x^, xj = min ^^^a, 

^'^' mi,bkCCk=x^-x, 

a^ > k=l, ..., 2 p. 



For any two fiindamental directions bp b/, e JR^ with 
j ^ k and j ^ k-\- p, let r(bp b/,) denote the cone in JR^ 
generated by bj and b/,, that is, r(bj, bi) = {x : x = bjaj-\- 
b,ai^aj,a,>0}. 

For any two points x^ and x^ e JR^, the vector x^ - x^ 
must be in some cone generated by two adjacent fiinda- 
mental directions, that is, x^-x^e r(bi,, b^+i) for some 
k= 1, ..., 2p. ThusjCo-JCj = 6jt^jt + ^jt+i^m for some non- 
negative scalars a/, and a^^+i. The vector x^ - x^ might 
also be in one or more cones generated by pairs of non- 
adjacent fiindamental vectors. However, the following 
Property shows that an optimal basis to the linear pro- 
gram LPD must correspond to adjacent fiindamental 
directions. 

Property 1: Suppose b^, bi, ^3 e E} are fiindamental 
directions of some unit ball B, Suppose x e T{bi, 62) 
with X = b^a^ + 62^25 and a^, o^ > 0. Suppose 63 e r(b^, 
62) with 63 = b,P, + ^2^. i3i + A > 1 and A, A > 0. If 
X e r(6i, 63) with X = b{Y^ + b^y^, y^, 73 > then Yi + 73 < 
a, + aj. If X e r(b^, b^ with x = b^y^ + bjyj, 73. 72 ^ 
then 73 + 72 < OTi + a2. 

Proof: Suppose x e r(bi, 63), then x = b^y^ + b^y^, = 
^i(7i + YsPi) + bjiYspj) by substitution. Since the repre- 
sentation of jc as a nonnegative linear combination of ^i 
and ^2 is unique, and x = b^a^ + 62^ then ai = 7i + ji^fiy 
and a^ = yA- Thus a, + a2 = r, + y^Pi + XsA = Xi + 
YsiPi + A) > 7i + 73- If-v e r(6i, 62) the proof is analo- 
gous, n 

Suppose that x^-x^e r(bf^ ^^t+i) for some k. Then the 
2 by 2 matrix [ b^ b^^^ ] is a feasible basis for the lin- 
ear program with 



^{h KiW^o-^d) 



and 



where e^ = (1, 1). 

Block distance may also be characterized in terms of 
the polar set B^ of the polytope B. The polar set B^ is 
also a symmetric polytope defined hy B^= {v \b\v<\, 
^ = 1, ..., 2/?}. In general, the facets ofB are in one-to- 
one correspondence with the extreme points ofB^. In 
R^, B^ has the same number of extreme points as B, 
which correspond to polar directions and are denoted 
by b\ for k= \, ..., 2/?. It may be shown that the polar 
directions are given by 
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Consider the dual of LPD, stated below as DLPD: 

DLPD: d^(x^,x^) =max v''(x^-xj 

s.t. blv < 1 k = \, ...,2p. 

The constraint set of DLPD is the polar set B^ which, by 
the Representation Theorem [5], may be written as B^ = 
{v = i:l^^X\^ ^l=A =1, ^ ^0}. Substituting into the 
dual objective function gives the equivalent characteri- 
zation of block distance in terms of polar directions: 

d^(x,,x,)= max bl^{x^-x,), 

k=l,...,2p 

Block distances are used to model travel distance in 
which the directions of travel are restricted to the fun- 
damental directions. The l^ distance is an example of a 
block distance with p = 2. Its fundamental directions 
are given by b^ = e^, ^2 = ^s ^3 = -^1? and ^4 = -^, where 
e^ is the ith unit vector in IR^. The polar directions for 
the /i distance are given by b^J = (1, 1), ^f = (-1, 1), 
^3 = -b], and bl = -b\, which are also the fundamental 
directions of the L distance, a block distance with p = 
2. 

Figure 1 illustrates a block distance with p = 4 and 
includes the fundamental directions bj,, k= 1, ..., 8, the 
unit ball and the polar directions b^ k = I, ..., 8. The 
components of bf, and ^^ used in Figure 1 are given 
belowiZ^f =(1,0XZ.[ =(f4), bl =(-hiX ^ =(-hjl 
Z.[ =(-1,0), bl =(-f,-fX Ih' =(^.-il li =(f. -f) and 

br=(m bi' =(i4x € -i-f.^'ii € -iMi €= 
(-1,-fxc =(-1 4x br =a-B^ r =0.-1). 

3. Linear Programming Formulations 

Ward and Wendell [8] presented two linear program- 
ming formulations of the one facility minmax location 
problem (with all Wi= 1) using block distance: one in 
terms of fundamental directions and one in terms of 
polar directions. These two formulations are given 
below for problem P2, and denoted as LPl and LP2. 
The expression of block distance in terms of polar 
directions is substituted into problem P2 to obtain the 
following: 

min z 

s.t. z > maxj^^ 2p ^Pf iPi-x) for / = 1, ..., 777. 

which is restated below as a linear program: 




bs^ bg 



bg 



Fig. 1. Fundamental directions by the unit ball, and polar directions 



LPl : min z 

s.t. z + wpf X > wpfpi for / = 1, ..., m, 
and k= \, ..., 2p. 

Substituting the expression of block distance in terms 
of fundamental directions into problem P2 gives the 
following: 



mm z 



s.t. z > iv.min^^^^a.^ for /=1,..., m 

Pi-^ = l^T=A^ik for/=l,...,m 



a,^ 



>0 



for /=1, ..., m, and k = l,...,2p 



which is restated as a linear program with all variables 
on the left as follows: 

LP2: min z 

s.t. ^-^i^/f^^(^ik -^ f^r /=1,..., m 

^^^li,^ic^iic=Pi for/=l,...,m 
a.^ >0 for / =l,...,m, 

and^ = l,...,2/> 
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The equivalence between the representations of block 
distance in terms of fundamental directions and in 
terms of polar directions implies that the formulations 
LP 1 and LP2 are equivalent. Only problem LP 1 will be 
considered in the subsequent development. 

The dual to LPl is given below where %, / = 1, ..., 
m, ^ = 1, ..., 2/?, are the dual variables. 



The direction 6? and the point p determine the 

r 1 1 

column , that enters the dual basis. The leaving 

b^A J 

column may be determined by using the simplex rules, 
that is by using the following equations to compute the 
components c/i, di and dz of the direction vector corre- 
sponding to the basic columns: 



r 1 1 




\ ^ 1 




\ ^ 1 




\ ~^ 1 


-,<_ 


d,+ 


_^A_ 


d,+ 


-,<_ 


d,= 


.-%*". 



Then the step size a and the leaving basic column are 
computed using the minimum ratio test: 



Since the dual constraints are of rank three, a dual basis 
has the form: 



\K, 



1 



1 



1 



w, b. w,- K n b. 



\ J 



and the dual basic variables are denoted by TZi^j,^, 71^^^/^, 
^/3,V where ij e {1, ..., m} and ^ e {1, ..., 2p} fory = 1, 
2, 3. The three weighted polar directions that determine 
a dual feasible basis are called basic weighted polar 
directions. 

Given a dual feasible basis, the dual simplex algo- 
rithm proceeds as follows. The basic weighted polar 
directions in the dual feasible basis correspond to active 
(equality) constraints in the primal, so that the variables 
z*, X* are determined by a solution to the following 
system of linear equations. 



"l 
1 




"z*" 

X* 


= 




1 


-.< J 




Wp, \ 



If z*, X* are primal feasible, that is, if 

z* >wpY(Pi-x'^), for all /= 1, ..., m, and^= 1, ..., 2p, 

then X* and z* are optimal. Otherwise, for some point 
p^ and some direction b% 

z*<wjbf{p,-x*) 

which implies that the point /?^ is outside the ball cen- 
tered at X* with radius z*/w^. Choosing the most violat- 
ed constraint corresponds to choosing a point of great- 
est weighted distance fromx*. 



a = mm < 



T'-"'- 



K, 



-dj* 



and the column 



>v,6" 



leaves the basis. The weight- 



ed polar directions in the new dual feasible basis are 
^/ *_ ^K_ 5 %^?5 ^i^d W/^ bl,^ . The algorithm continues 
until primal feasibility is achieved. Problem LPl is 
bounded and feasible, so that an optimal solution exists. 



4. Geometric Interpretation of tlie Dual 
Basis Update 

The update of the dual basis in the simplex algorithm 
applied to problem LP2 is analyzed in terms of the 
geometry associated with the basic weighted polar 
directions w,i>?. 

Ij Kj 

Consider a dual feasible basis 



1 



1 



W; b, w. b. 



1 



with basic dual variables n^^^j,^, 71^^^/^, n^^^k^, for ij e {1, ..., 
m}, kj e {1, ..., 2p} andy" = 1, 2, 3. Assume the weight- 
ed polar directions in the dual feasible basis are ordered 
counterclockwise with respect toy = 1, 2, 3. We adopt 
the notation that ify = 1 , theny -1=3, and ify = 3, then 
7+1 = 1. 

Geometrically, a dual feasible basis implies that the 
vector may be expressed as a convex combination of 
basic weighted polar directions, that is. 
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"%b? 



'^i,,h. 



If the dual feasible basis is non-degenerate, then 
71, J, > fory = 1, 2, 3 and is a strict convex combina- 
tion of the basic weighted polar directions. In this case 
each basic weighted polar direction is contained in the 
cone generated by the negative of the other two basic 
weighted polar directions, that is, Wi,bl, e 
int(r(-wM, -W;blJ) for eachy = 1, 2, 3, where 



int{T{-w,b\,-w,^y,J) = 



{x: x = -W;bl.a- 



bl 



w,,bi,fi\a,p>0}. 



Also, for any weighted polar direction wjbl, there is 
some / = 1, 2, 3 so that wj)^ e Ti-Wibl, -w-. bl ). 
Figure 2 illustrates the non-degenerate case with solid 
arrows corresponding to the basic weighted polar direc- 
tions Wi,bl, and dashed arrows corresponding to ->V//»L 
fory" = 1, 2, 3. For the non-degenerate case the basic 
weighted polar directions form a simplex in R^. Figure 



2 also illustrates w^b^ e r(-Wi,bl, 



-Wi bl ), 

V+i Vi 



% .l^J.L ^ 



"''v+'bg^^, 




>^^,K, 



Fig. 2. Non-degenerate basic feasible weighted polar directions. 

In the degenerate case, 7i, jt = for exactly oney = 1, 
2, or 3, and %-. t > 0, 7i, ^ > 0. Figure 3 illustrates 
the degenerate case. Note that if 71/ ;;, = 0, then Wi,^bl, ^ e 
^{-^ibl. -^i bl ), and w, bl e' r(->v,6?, -w^ bf ), 
but Tv,.^^ ^ r(-Wi bl , -TV,- 6^ . Also, if 71,- ^ = 0, then 
for any weighted polar direction w^bl, either w^bl e 
T(r^M, -w:bl) or >v,6? e n-wM, -w,bl% or 
>vJ?G R'\ (n-wM, -w,blJ<jTi-w,bl, -w,bl)}. 






-^i.K, 




^iX, 



Fig. 3. Degenerate basic feasible weighted polar directions. 



y, and since Wi,_bl^^ e ^-^ifil., -^ij,bl.J, then tv^^? 
replaces tV/ 6^ . The new basic weighted polar direc- 
tions are tv^6? , w^bl, and w^^^bl,^^. If w^^? e mt(T(-Wibl,, 
-TV, ^? )), then the new basic weighted polar directions 
are non-degenerate. If w^b^ coincides with -Wibl, or 
with -TV/ bl, , then there is a tie for the replaced 
weighted polar direction and the new basis is degener- 



ate. That is, if w^b^ coincides with -w^bl,, then w^bl e 
r(-w,bl,-w, bl )mdwy,G r(->v,4/->v, bl ) and 
either tV/ 6^ or tV/ 6^ may be replaced. If wjb^ coin- 



cides with ->v, bl , then>v„^?e r(->V:^?,->v, bl ) and 
wjbl e r(-TV,- 6? , -TV; ^2 ) and either w-. bl or tv,-62 
may be replaced. 

For the degenerate case suppose 71, jt = 0. If tv^6? e 
r(->V/,^?,, -^/,+i^Vi)' then >v^6?replaces^ ^6?^^. The new 
basic weighted polar directions are w^, w^bl, and 
w^^b^^^ Ifw^, G int{T{-w,bl,, -w^jl.J% then the new 
basic weighted polar directions are non-degenerate. If 
w^ coincides with -w^bl,, then tv^^J e r(-TV,/>^, 
-TV, bl ) and tv„6? e r(-Wibl, -tv, 6? ) and either 
TV, 6^ or TV, ^t may be replaced and the new basis is 
degenerate. However, if wjb^ coincides with -tv,^^^^^^, 
then Wqbl e T(-Wi,bl,, -Wi^pl,^^ only, so that only 
TV, bl_ is replaced. The new basis is degenerate. 

The degenerate case with 71, ;;, = and wjbl e 
r(-TV,/>^ , -TV/ bl, ^), is analogous to the degenerate case 
in the preceding paragraph withy + 1 interchanged with 
y - 1 throughout. 

Finally, consider the degenerate case with n^ /, = 
but w,b', € T{-w,fy\, -^lAJ, and w,bl € T(-w,bl.. 
~^i_fik._)- In this case tv^^^ replaces w^bl. The new basic 
weighted polar directions are tv,^^^? ^, tv^^? and ^i.^fil.^^ 
and remains degenerate. 

The following Property shows that the new weighted 
polar directions determined by the geometric proce- 
dures above are basic feasible. 



For the non-degenerate case, the basis update rule is Property 2: Suppose the weighted polar directions 



given as follows: if TV^^^e r(-TV,/>^, -tv, bl, ), for some 



TV, bl , wM , and tv, bl , are basic feasible. For either 
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the non-degenerate or degenerate case, suppose tv^^? = 
-w^fylfi - w,.Jl^ Pj^, with Pj, Pj,, > 0. Then the weight- 
ed polar directions tV/ bl, , w^b^, and Wibl, are basic fea- 
sible. 

Proof: The assumption w^b^ = -WiPl^j - tV/ b^. Pj+i 
with Pj, Pj,, > implies wMjPj + wA + ^JVP^ = 
O.Letr = j3,+ l+j^,i>O.Let7r;, = ^>0,<,= ^> 
0, K^^i,^, = ^ > 0. Thus tt;,. + tt;, + tt;,! ..^1 = l, and 
the weighted polar directions tv^ bl_ , w^b^, and Wibl_ 
are dual feasible. Observe that if either j3^ = or Pj+]^ = 
0, the new basic weighted polar directions are degener- 
ate, n 



Dividing through the last equation by the right hand 
side, and comparing the resulting two equations to (1) 
and (2), shows that 



?r,. 



-d,., 



''-''''-' A. +A . +1 



SO that 4_ < since 7t, jt ^ 0. Equation (1) implies 
that 



-h;. bl =w.bl 



,0 J'J -L-.., aO 'y+i'Vi 



71. 



i t 






?r, 



'i Ir 



The next property shows that the geometric replace- 
ment rule corresponds to the minimum ratio rule of the 
simplex algorithm. 



Property 3: Suppose the weighted polar directions 

Wi bl,_, WiPl,, and tV/ bl, , are basic feasible, and sup- 
pose wj)l is the entering weighted polar direction and 
that Tv„^5 e r(-Tv.^? , -TV, bl ). The geometric rule that 
Wqb^ replaces tv^ bl,_ is equivalent to the minimum ratio 
rule of the simplex algorithm. 

Proof: The proof is given for the non-degenerate 
case. The proof of the degenerate case is similar. First 
we show that if wM e T{-Wibl, -w^ bl ), then 

Dual feasibility implies the following equations: 



^i bl %; k + 1^/^?^/ Jt + 1^/ ^? ^/ ^ =0 (1) 

^/ k ^ T^i k^ ^/ jt = 1- (2) 

The components of the direction vector determined by 
the vector w^bl are given by the following equations: 



«>r 



dj- 



4.1= !■ (4) 



The assumption that w^b^ e r(-Tv,^^ , ->V/^^^?^^) implies 



-■^iblPik-^i bl A- >t =>vX 

7 J f J 7+1 7+1 J+l'Vl ^ 



(5) 



where j3, j, > and Pi.^^^t^^ > 0. Equations (3) and (5) 
combine to give 

and equation (4) implies that 

-4-1 + (A,,, vi - 4) + (^..,^,.1 - 4+i) = A.., + /3^.,,.,,, + 1- 



Substitute this expression into equation (3), multiply 
through by %i k. and divide through by -^-i to get 



n. 



'j^^j-^ A \... aO 



n, 



'vi'^^1 



{K, , + ' ' d^w.bl +(;r, , + ''' '' d^^,) 

"y-i "y-i 



)>V Z'' 



;r, 






(6) 



If c^ > 0, the first coefficient is positive. If dj < 0, the 
first coefficient is non-negative if and only if 



^i , k , ^i k 



-dj_, -dj 



A similar argument holds for the second coefficient. 
The third coefficient is positive. Thus, the new coeffi- 
cients are non-negative if and only if a = '^"^'''"' is the 
minimum ratio. 

To prove the converse, suppose that min^i 23(^7^ '- 
4 < 0} = ^^, and show that w^bl e T(-wM. 

-w, bl ). ^"' 

Vi j+y 

By the minimum ratio assumption, dj_i<0. Ifdj < 0, 
then 'CT'~' < -zj- by the minimum ratio, and the first 
coefficient in equation (6) is non-negative. Ifdj > 0, the 
first coefficient in equation (6) is non-negative. A sim- 
ilar argument holds for dj^i and the second coefficient in 
equation (6). Thus wfi^, e T{-W;bl,, -Wi,^bl,J. n 

Properties 2 and 3 provide a geometric rule that 
could be used to determine the leaving column of a dual 
basis in the basis update step. However, the equivalent 
minimum ratio rule of the simplex algorithm is more 
efficient. 

In the next section, the equivalence is used to show 
that the minimum ratio rule may be used to update the 
dual algorithm applied to the Euclidean distance one- 
center problem, which is an improvement over existing 
geometrical update rules. 
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The results of the last two sections may be extended 
to polyhedral norms with consideration given to the 
asymmetry. The preceding results may also be extend- 
ed to ^" with additional complexities and notation. 



5. Euclidean Distance Problem 

Considered next is the one-center location problem 
using Euclidean distance and assuming all weights Wi = 
1 . The constrained version of this problem is written as 
follows: 

P2: min z 

sJ. z > l2(Pi, x) i = 1, ..., m. 

Euclidean distance may be considered as a block dis- 
tance with p —> oo and whose unit ball is the circle of 
radius 1. There is a fundamental direction for every 
point on the unit circle and each fundamental direction 
is of unit length. The polar directions for Euclidean dis- 
tance are identical to the fundamental directions. Given 
two points Pi and x, the unit vector ^^^ represents the 
polar direction from x to /?,. The polar direction that 
maximizes the dot product with/;, - jc is the unit vector 
^z^. Thus the Euclidean distance may be interpreted in 
terms of polar directions and is written as: /sQ^,, x) = 
-^Tt\(Pi ~ ^)' Problem P2 is then written as: 

P2: min z 

An alternative and equivalent statement of problem 
P2 is to use the squared Euclidean distance, that is, l^ip^, 
x) = ipi - xY ipi - x). The squared Euclidean distance 
between the points Pi and x may also be considered as 
a block distance with/? -^ ^ and whose unit ball is the 
circle of radius \Pi-x\. There is a fundamental direction 
for every point on the unit circle, and the polar direc- 
tions are identical to the fundamental directions. Given 
two points Pi and x, p^ - jc is the polar direction from x 
topi. Problem P2 for the squared Euclidean distance is 
written as: 

P2: min z 

sJ. z>{pi-xY{pi-x) /=1, ..., m. 

The squared Euclidean distance problem has the advan- 
tage that the constraints are differentiable. This leads to 
the Karush-Kuhn-Tucker conditions for problem P2: A 
point X and a radius z are optimal to P2 if and only if 
there exists tt, for / = 1, ..., n, such that: 



z > {p,-x) (p,-x) 


i=l, .. 


., m 


(V) 


I>.=1 






(8) 


X^.U-^K-o 






(9) 


;r,. >0 


/=!,.. 


, m 


(10) 



(z-(/,,-xf(/,,-x));r,=0 i=l. 



m 



(11) 



Conditions (7) represent primal feasibility, condi- 
tions (8), (9), (10) represent dual feasibility, and condi- 
tions (11) represent complementary slackness. The dual 
feasibility conditions are analogous to constraints of the 
dual problem of LPl. For a given jc, the dual constraints 
are of rank three and a dual basis has the form: 



1 1 1 

-X p,-x p,- 



The dual variables are denoted by %i , tt, and %i respec- 
tively. The polar directions Pi - jc for y = 1 , 2, 3 that 
determine a dual feasible basis are called basic polar 
directions. 

Elzinga and Heam [3], report an efficient dual based 
algorithm for solving the Euclidean distance min-max 
problem that is described as follows. Choose three 
points/;,, ij e {1, ..., m} andy = 1, 2, 3. In the primal 
phase of the algorithm, the min-max solution x for 
these three points is determined as follows. If the three 
points form an acute triangle, then x is the intersection 
of the perpendicular bisectors of the line segment 
between any two pairs of points. If the three points 
form an obtuse triangle, then x is the mid point of the 
line segment between the pair of points forming the 
longest side (opposite the obtuse angle) of the triangle. 
In either case, z is the maximum distance from x to the 
three points. 

Next the solution x and z is checked for primal feasi- 
bility. If the distance between x and Pi does not exceed 
z for all / = 1 , ..., m, then x and z is an optimal solution. 
Otherwise, some point p^ is chosen that violates primal 
feasibility. The Elzinga-Heam algorithm provides a 
geometric procedure to determine which of the three 
points/;, ,7" = 1, 2, 3 is replaced hy p^ and the algorithm 
continues. 

The Elzinga-Heam geometric procedure for deter- 
mining which point is replaced by/;^is described as fol- 
lows: Determine the point/;, ,7" = 1, 2, 3 that is of great- 
est distance from/;^ say/;,. If the point/;, ^ (or/;, ^) and 
/;^are on the same side of the line through/;, andjc, then 
/;^ replaces /;, ^ (or/;, ^). Figure 4 illustrates an example 
where p^ replaces /;, . 
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Fig. 4. Euclidean distance example. 

In terms of dual feasibility, the three points Pi fory = 
1, 2, 3 correspond to the basic polar directions/?,- -jc in 
the dual basis. If the three points form an acute triangle, 
the corresponding basic polar directions are non-degen- 
erate and form a simplex in IE?. If the three points form 
an obtuse triangle, the corresponding basic polar direc- 
tions are degenerate. The point p^ corresponds to the 
polar direction/?^- jc that will become basic and replace 
some existing basic polar direction. Figure 4 illustrates 
the basic polar directions. 

The Elzinga-Heam procedure is now compared to 
the geometric replacement rule developed above. 
Suppose p^-x e ^(-Pi + -V, -/?,-^^ + x). Then the point 
farthest from/?^ is either/?,- *^^/^^+^ ^y ^^^^ feasibility 
/?,- ^ - jc e r(-/?,- + X, -pi,^^ + x). The replacement rule 
developed above states that/?^-Jc replaces/?,-^ -x. The 
inclusion/?,-^ -jc e r(-/?,- + jc, -p, ^ + jc) imphes that/?,- ^ 
and/?^ are on the same side of the line through/?,- and x, 
and that /?,-^ and /?^ are on the same side of the line 
through /?,;^ and x. Thus, the replacement rules are 
equivalent. If the basic polar directions are degenerate, 
a similar argument shows the equivalence of the 
replacement rules. 

Since the geometric replacement rule for block dis- 
tance is equivalent to the minimum ratio rule of the 
basis update, the minimum ratio rule provides an alter- 
native to the Elzinga-Heam geometric rule for selecting 
the point to be replaced in the dual based algorithm for 
the Euclidean distance min-max location problem. 

The algebraic replacement rule for the un-weighted 
Euclidean distance min-max problem has been extend- 
ed to the weighted Euclidean distance min-max prob- 
lem in a subsequent paper by the first author. 
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